Full optical responses of one-dimensional metallic photonic crystal slabs 
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We reveal all the linear optical responses, reflection, transmission, and, diffraction, of typical 
one-dimensional metallic photonic crystal slabs (MPhCS) with the periodicity of a half micrometer. 
Maxwell equations for the structure of deep grooves are solved numerically with good precision 
by using the formalism of scattering matrix, without assuming perfect conductivity. We verify 
characteristic optical properties such as nearly perfect transmission and reflection. Moreover, we 
present large reflective diffraction and show that, in the energy range where diffraction channels are 
open, the photonic states in the MPhCS originate from surface plasmon polaritons. 

PACS numbers: 42.70.Qs, 42. 25. Fx, 42.25.Bs, 73.20.Mf 



Electromagnetic (EM) waves in artificial structures 
were one of the old problems^ and had been thought to 
be already resolved. Since the observation of extraordi- 
nary transmission of perforated silver film in 1998, 3 opti- 
cal responses of metallic photonic crystal slabs (MPhCS) 
have attracted much interest as a novel phenomenon in- 
volving surface plasmons. Many researchers have been 
stimulated by the phenomenon and have tried to explain 
the properties of transmission in MPhCS. 

Even when we restrict an object to a one-dimensional 
(ID) MPhCS as drawn in Fig.^ the zeroth-order trans- 
mission was already discussed in not a few reports 4 *^ 
and explained qualitatively. In order to evaluate 
the transmission numerically, Maxwell equations were 
solved with additional assumptions such as perfect 
conductivity*^ In those reports, transmission coefficient 
was focused on and calculated. 

On the other hand, diffraction has not been fully ana- 
lyzed, to our best knowledge, although the surface plas- 
mons in MPhCS are usually connected to diffraction. In 
fact, diffraction in optically thick metallic grooves such as 
Fig- H has been a tough subject since the first era of sur- 
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FIG. 1: Structure of ID MPhCS; the periodicity is 500 nm, 
the width of air slit is 100 nm, and the thickness d is 300 
nm. The MPhCS is on the substrate of quartz. Incident light 
(arrows) travels in the xz plane. Incident angle 9 is the angle 
between the normal and incident directions. The incident 
light is p-polarized. 



face plasmons in 1980s. Maxwell equations for metallic 
grooves were analytically solved under the assumption of 
perfect conductivity;- however, the results based on the 
hypothesis deviate from experimental results quantita- 
tively. Thus, there has been strong motivation for de- 
scribing the optical responses of MPhCS quantitatively. 

The formalism to calculate full optical response 
including diffraction was tried to construct, inde- 
pendently of the extensive studies on extraordinary 
transmissioni&iLiSiiiiiSiiii To include diffractive com- 
ponents, scattering-matrix (S-Matrix) formalism is 
suitable!?'?' 1 ?: 11 ' 14 The S-Matrix method is a kind of 
transfer-matrix method for the propagation of EM wave. 
The method was invented to ensure numerical stability?^ 
We adopt the S-Matrix formalism; it is based on the ordi- 
nary local response of EM fields, and the explicit expres- 
sion for optical responses was presented in Ref. Il4l The 
formalism worked quite well for thin ID gold PhCS^ 

In this Communication, we solve Maxwell equations 
for the ID MPhCS accurately and evaluate all the lin- 
ear optical responses. In particular, we explore the large 
diffraction and seek the physical insight. Furthermore, 
we discuss the physical meaning of this computation. 

Figure n shows the structure we mainly examine here. 
Rectangular silver rods are infinitely long and periodi- 
cally located at the 100-nm interval on the quartz sub- 
strate. Incident light is p-polarized, that is, the polar- 
ization of electric field is perpendicular to the axes of 
Ag rods. Setting x axis to be the periodic direction, the 
electric field of p-polarized incident light is expanded as 



E(r,i) = E m exp(ik m ■ r - iwt), 

m 

k x ,m = k x + 27rm/a, (m = 0, ±1, ±2 
k x = &osin#, 
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where a is the periodicity, fco is wavenumber of incident 
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FIG. 2: (a) Reflectance and (b) transmittance spectra of the 
ID MPhCS at incident angle = 0°. 

light, 8 is defined in Fig. ^ and k y = 0. The electric 
fields in the MPhCS and substrate are expressed accord- 
ingly. The magnetic fields are described likewise. In the 
present S-Matrix formalism, the infinite expansion of Eq. 

gives the Fourier-based linear equations which are ex- 
actly equivalent to Maxwell equations. Practically, the 
truncation of expansion is inevitable. We cut Eq. (|TJ 
at 77i = ±TO max and then the total number of harmon- 
ics -/Vtotai is 277i max + 1. Numerical computations were 
executed in the highly vectorized supercomputers, which 
enable to make 32 parallelization. The complex dielectric 
constant of silver is taken from Ref. 

Figure [21 displays (a) reflectance and (b) transmit- 
tance spectra calculated at incident angle = 0° and 
-^totai — 1265. Diffraction is calculated simultaneously 
in this method. The fluctuation of the calculated values 
is within 1% for AT total > 1200 in this energy range. 

Nearly perfect transmission at 1.2 eV in Fig. |2{b) is 
due to the waveguide mode in the air slit. The qualitative 
analysis of the mode was reported previously^*£*& and the 
Fabry-Perot nature was pointed out. The property is 
reproduced in our calculation. The waveguide mode does 
not appear above 1.7 eV where diffraction is effective, 
and transmission is not strongly enhanced. Reflectance is 
nearly 100% at 1.63 and 2.36 eV in Fig.^a); the response 
is called Wood anomaly^ It is usually observed at the 
resonances of surface plasmons in metallic grooves. 

We note here the reason why the large iVtotal is re- 
quired. This is mainly because there exist the resonances 
under p polarization, due to the waveguide mode and 
surface plasmons. In a good metal such as silver, the res- 
onances induce the spatially modulated and intense local 
EM fields in and near the MPhCS; the situation is in gen- 
eral hard to express with a small ./Vtotai- Therefore, reli- 
able numerical calculation for optical resonances requires 
the large -/Vtotai and actually becomes quite demanding. 
However, we can achieve optical responses with good pre- 
cision by using a large -/Vtotai in a wide frequency range 



FIG. 3: (a) Solid line: total diffraction spectrum of the ID 
MPhCS. (b) Dashed line: transmissive diffraction. Dotted 
line: reflective diffraction. The sum of these diffractions is 
the total diffraction. The incident angle 9 is 0° . 



from near infrared to ultraviolet. On the contrary, in the 
off-resonant frequency range, far less -/V to tai is enough to 
obtain the solution of Maxwell equations. Indeed, we 
can calculate the optical responses within 0.2% fluctua- 
tion by using a -/V to tai of a few tens under s polarization, 
when the polarization of electric field is parallel to the 
rod axes. 

Figure 02a) presents total diffraction spectrum under 
the normal incidence (solid line). We use the conven- 
tional notation in diffraction theory; R n means the nth 
reflective diffraction and T n the nth transmissive diffrac- 
tion. Then, reflectance is Rq, transmittance is To, and 
the total diffraction is Xm^o(-^n + T n ). Dashed line in 
Fig. |2Ib) shows transmissive diffraction ^2 n ^ T n . Sim- 
ilarly, the reflective diffraction J2 n ^o ^ n ^ s represented 
with dotted line. 

The components of diffraction have the following prop- 
erties. The transmissive diffraction T±\ emerges at 1.70 
eV. The sharp onset corresponds to the opening of the 
first-order channel of Bragg diffraction at the interface 
between the quartz and MPhCS. The second channel 
opens at 3.40 eV and induces T±2- On the other hand, 
the reflective diffraction appears above 2.48 eV and is 
composed of R±i in Fig. EJb). The onset is due to the 
opening of the first-order channel of Bragg diffraction at 
the interface between the MPhCS and air. The sum of 
R±i reaches the maximum of 70% at 2.75 eV. The spec- 
tra of diffraction are discussed later. 

To examine the dispersion of resonances above 1.3 eV, 
we show extinction spectra at incident angle 9 = to 50 
degrees in Fig.^Ja). Extinction is — In (To). The peaks of 
the extinction spectra are the minima of transmittance 
by definition and plotted with solid circles in Fig.^Jb). 

We compare the dispersion of extinction with that of 
surface plasmon polariton (SPP). SPP at flat interface is 
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FIG. 4: (a) Extinction spectra of the ID MPhCS at inci- 
dent angles from to 50 degrees, (b) Dispersion relation 
of extinction peaks (solid circles). Dotted and dashed lines: 
the dispersions of reduced SPPs at the air-MPhCS and the 
quartz-MPhCS interfaces, respectively. The periodicity a is 
500 nm. The wavenumber k x is described in the text after 
Eq. ©. 



generally expressed by 

where w is the frequency of incident light, c is the ve- 
locity of light in vacuum, and fcii is wavenumber of SPP 
along the x axis, e is dielectric constant of air or quartz, 
and £Ag is the complex dielectric constant of silver. Note 
that when £Ag is complex number, fcii is also complex. 
For silver in the energy range of interest, Eq. @ holds 
with replacing e^g by Re(eAg) and k\\ by Re(fc||)iiS We 
set k x = Re(fc||) in accordance with Eq. QJ. The dis- 
persion of SPP is evaluated and shown in Fig. 0fb) after 
reducing into the Brillouin zone. The dotted and dashed 
lines represent the SPPs at the air-Ag and the quartz- Ag 
interfaces, respectively. 

The dispersion of extinction is in good agreement with 
the dispersion of reduced SPPs. The agreement indi- 
cates that the reduced SPPs are resonantly excited, and 
strongly suggests that the photonic states in the MPhCS 
have formed from the reduced SPPs above 1.3 eV. The 
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FIG. 5: Total diffraction spectra at various thickness d. The 
periodicity and the width of air slit are shown in Fig. The 
incident angle 9 is 0°. Arrows indicate the opening of nth 
diffraction (n = ±1,±2) and the resonances of SPPs. For 
example, SPP (1Q) and SPP (1A) stand for the first-order res- 
onances at the quartz-MPhCS and air-MPhCS interfaces, re- 
spectively. The energies are taken from Fig. |ljb). The spectra 
at d > 50 nm are offset for clarity. 



agreement is also consistent with the previous analysis 
with the method of rigorously coupled wave^S 

As shown in Fig. the sharp onset of diffraction at 
1.70 eV corresponds to the opening of T±\. On the other 
hand, as shown in Figs.|21and0J the sum of R±i gradually 
increases and has the maximum at 2.75 eV; the maximum 
is close to the resonance of the second-order SPP at 2.80 
eV and k x = 0. How can we explain the enhancement 
of diffraction? To extract more information, we examine 
the diffraction with varying only the thickness of MPhCS. 

Figure represents total diffraction spectra at the 
thickness d = 20, 50, 100, 300, 500, and 800 nm under the 
normal incidence. Arrows indicate the energy positions 
at which the channels of diffraction open, and the reso- 
nances of reduced SPPs. SPP( IlQ > and SPP( nA > denote 
the nth resonances at the quartz-MPhCS and air-MPhCS 
interfaces, respectively. The positions of SPP resonances 
are taken from the dispersion in Fig. 0Jb). In all the 
spectra, diffraction arises at 1.7 eV and is comprised of 
T±i. At 1.7-2.4 eV, diffraction reaches a few tens of per- 
cents and the spectra show similar tendency irrespective 
of the thickness. Above 2.4 eV, the diffraction spectra 
depend on the thickness. In the thin MPhCS of d = 20 
and 50 nm, the components T±\ and R±i coexist, and are 
not prominently enhanced at 2.7-2.8 eV. In the optically 
thick MPhCS of d > 100 nm, large peaks of diffraction 
appear at 2.7-2.8 eV and are dominantly composed of 
R±i; see Fig.Efb) for the case of d = 300 nm. The peak 
positions change with the thickness, while the extinction 
spectra show that the energies of SPP resonances are 
common in the MPhCS of d > 100 nm. Thus, though 
the large peaks at 2.7-2.8 eV are close to the resonance 
gpp(2Q) a t 2.80 eV, they cannot be simply ascribed to 
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the resonant enhancement by reduced SPPs. 

By varying the thickness from d = 20 to 800 nm, one 
can observe that the diffraction spectra change in a quite 
complicated way. The changes in the diffraction spectra 
probably come from the changes of the photonic states 
in the MPhCS. We believe that the prominent peaks of 
diffraction indicate the resonant emission modes in the 
MPhCS. 

To find the resonant modes in the MPhCS, it is a key to 
solve the equation of det(<S) = 0, where S is a S-Matrix 
for the ID MPhCS, satisfying \in) = S\out), and \in) 
and \out) are input and output states, respectively, ft is 
surely demanding and a future issue to solve det(<S) = 
for a large iVtotal- The solutions of the equation and the 
photonic modes were obtained for the two-dimensional 
PhCS made of dielectrics with a small iVtotai^^i 

We have obtained a numerically accurate solution of 
Maxwell equations for the ID MPhCS and have shown 
the optical responses. We would like to compare the 
present results with measured ones quantitatively. For 
the precise comparison of the present calculated results 
with measurement, it is better to compute with using the 
complex dielectric constants obtained from the material 
used in the measurement, because the values can change 
under the condition of fabrication and environment. 22 
Such comparison is a test for the validity of the ordinary 
local response of EM fields such as the electric flux den- 
sity D(r) = e(r)E(r). If there exists definite discrepancy 
beyond the errors in the numerical calculation and ex- 
periment, it would be necessary to take into account the 
microscopic nonlocal response of EM fields^ It is an in- 
teresting subject to find the metallic nano-structures in 
which the nonlocal responses emerge distinctively. 

In the S-Matrix formalism, one can obtain absorption 
A by the relation of A = 1 - ^ aU „0Rn + T n ). Thus, 
the absorption can be evaluated from all the compo- 
nents in optical measurement. The absorption spectrum 
is also measured directly by photoacoustic methodjSi 2 ^ 
the method is useful for diffractive media. 

At the end of discussion we mention the previous re- 



ports. Nearly perfect transmission, which appears at 1.2 
eV in Fig.|2Jb), was attributed to the waveguide mode in 
the slit mainly from the distribution of EM fields^ The 
high-efficient transmission was also analyzed under the 
assumption of the perfect conductivity only in the slit*£ 
Besides, the coupling of SPP with light was recently ana- 
lyzed with extracting a effective term of S-Matrix and the 
transmission was explained qualitatively. 2 ™ In this Com- 
munication, Maxwell equations for the ID MPhCS have 
been solved numerically without any further technical as- 
sumption except for the inevitable truncation. From the 
solution, the optical responses have been evaluated di- 
rectly in a wide frequency range from infrared to ultra- 
violet. 

In conclusion, we selected typical ID MPhCS made of 
silver and have computed the full optical responses pre- 
cisely using the measured complex dielectric constant in 
the S-Matrix formalism. This is the first achievement 
concerning realistic calculation for ID MPhCS with the 
structure of deep groove. We have confirmed the res- 
onant excitation of reduced SPPs and found that the 
large reflective diffraction appears in the optically thick 
MPhCS. The origin has been discussed in view of res- 
onant enhancement by reduced SPPs and of resonant 
emission modes in the MPhCS. The present computa- 
tion, assuming the local response of EM fields, can be 
readily compared with measurement. 
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